Atom-molecule equilibration in a degenerate Fermi gas with resonant interactions 
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We present a nonequilibrium kinetic theory describing atom-molecule population dynamics in a 
two-component Fermi gas with a Feshbach resonance. Key collision integrals emerge that govern the 
relaxation of the atom-molecule mixture to chemical and thermal equilibrium. Our focus is on the 
pseudogap regime where molecules form above the superfluid transition temperature. In this regime, 
we formulate a simple model for the atom-molecule population dynamics. The model predicts the 
saturation of molecule formation that has been observed in recent experiments, and indicates that 
a dramatic enhancement of the atom-molecule conversion efficiency occurs at low temperatures. 



Over the past year, extraordinary progress has been 
made using a Feshbach resonance to control the atomic 
interactions in trapped Fermi gases By adjusting 

an external magnetic field, the energy of a closed-channel 
bound state e TCS (B) can be tuned relative to the k = 
threshold of a colliding pair of atoms Q. This has 
two dramatic effects on the system: atoms interact very 
strongly and molecules can be produced. 

Due to the resonant state coupling, the atom-atom 
scattering length a(B) can be tuned to very large val- 
ues, either positive or negative. One of the most stun- 
ning applications of this remarkable control of atomic 
interactions has been the recent achievement ofpairing 
superfluidity in a dilute gas of fermionic atoms |2j. This 
remarkable tool allows one to continuously traverse the 
so-called BCS-BEC crossover region |4j. 

In this paper we develop a kinetic theory for a reso- 
nantly coupled atom-molecule mixture in the pseudogap 
regime relevant to current experiments pj. A key in- 
gredient of our theory is a pair of newly identified collision 
integrals describing how the atom and molecule distribu- 
tions come into chemical and thermal equilibrium. As a 
first application, we derive a non-interacting model for 
the population dynamics that predicts the saturation of 
molecule formation observed in recent experiments 0,0- 
It also predicts that the efficiency of atom-to-molecule 
transfer increases as the temperature is lowered. 

During resonant collisions pairs of atoms form tran- 
sient molecules, without the influence of the third atom 
required for recombination in nonresonant collisions 
Such molecules are typically in a high-lying vibrational 
state that is very susceptible to collision-induced decay 
to a lower bound state. The energy released in such colli- 
sions quickly heats up the g limiting factor that has 
frustrated attempts to achieve Bose-Einstein condensa- 
tion of Feshbach molecules made from bosonic atoms . 

For fermionic atoms such inelastic collisions are sup- 
pressed by approximately three orders of magnitude due 
to Fermi statistics Q, which results in lifetimes of the 
gas sample long enough that the molecules may relax 
to equilibrium in the trap. Consequently, an important 
question is how to properly describe the atom-molecule 
population dynamics in this regime 0, Il0| . Theories 



based on a two-atom description, from which a two-level 
Landau-Zener treatment has emerged as a useful model 
should only be valid in the regime t sc « ( < t co \, 
where r sc ~ a/v is the duration of a scattering event and 
t co \ ~ 1 /nav is the average time between collisions. Here, 
a is the s-wave scattering length, v is the mean velocity 
of an atom in the gas, and a is the collision cross-section. 

In order to describe the current experimental regime, 
Tcoi <C t, we have derived coupled kinetic equations for 
the atom and molecule distributions above the superfluid 
transition temperature using the Kcldysh nonequilibrium 
Green's function formalism (T^El • Based on a self con- 
sistent diagramatic expansion, the atom and molecule 
self energies are treated on the same footing, leading to 
a renormalization of the atom and molecule mean field 
energy shifts. Before discussing the salient features of 
the theory, we first give a brief summary of the kinetic 
equations; a detailed derivation will appear in a future 
work 14] . We then derive a simplified model describing 
the atom-molecule population dynamics, which we use to 
understand recent experiments. 

The starting point of our theory is the coupled boson- 
fermion Hamiltonian H 



H = 



dr 



V^(r)#4r)^(r) + ft (r)H M (r)ftr) 



+ n dr 



ft (r)^ (r)fr (r) + Tp\ (r)$J (r)ftr) 



g A / drV'j(r)V'{(r)'0|(r)^ T (r), 



(1) 



where a sum over spin states a = {Til} is implied in the 
first term. The atom ^(r) and molecule ftr) field opera- 
tors obey Fermi and Bose commutation relations, respec- 
tively. The single-particle hamiltonians for the atoms and 
molecules are given by Ha(t) — — (S 2 /2m)V 2 -I- Ua{?) 
and H M {r) = -{h 2 /4:m)X7 2 + U M {r) + e rcs , where U A (r) 
and U m (r) are the external trapping potentials for atoms 
and molecules and m is the atomic mass. In experiments, 
the energy e Tes (B) of the resonant molecular state can be 
tuned by adjusting an external magnetic field B. The 
coupling to the resonant bound state is described by the 
two terms in the second line of with coupling constant 



2 



k; the purely open-channel part of atomic collisions is de- 
scribed by the last term with the interaction constant g&- 
Starting with the many-body Hamiltonian QJ, we de- 
rive equations of motion for the atomic and molecular 
distribution functions / CT (p,r,i) = J dxe tp ' x/ ' ft {'0j.(r + 
x/2)^(r - x/2)) t and f M (p,T,t) = J dxe^/^r + 
x/2)i/)(r — x/2))t, where (• • -) t indicates a nonequilibrium 
statistical average. The steps leading to the kinetic equa- 
tions are rather involved and are beyond the scope of this 
letter 0,0,0; here we state the results. The kinetic 
equations describing the atom-molecule dynamics are 



dt 

dfM(p,r,t) 
dt 



{e (T (p,r,t),f a (p,r,t)}=l„+I a ,{2) 
{e M (p,r,t), f M (p,r,t)} =2 M , (3) 



where {A, B} = V r A ■ V p B - V p A ■ V r B is 
a Poisson bracket. The time-dependent renor- 
malized energies for the atoms and molecules are 
given by (with the position dependence suppressed) 

e a (p) = e^(p) + Pfdp'g A Ap,p')f-*(p')/h 3 + 
V J dp'g K ,<r(p,p' - p)/a/(p')/^ 3 and e M (p) = £m (p) + 
Vf dp'g^ip - p', p')[/ T (P - P') + /|(P') - 1], where V 
indicates the principal value. The noninteracting parts 

of the energies are given by e^?(p, r) = p 2 /2m + Ua{y) 
and e^(p,r) = p 2 /4m + Um(t) + e res . The renor- 
malized atom-atom interaction strength is <?a,o-(p, p') = 
9A+9nAv,P') where g K ,a{p,p') = K 2 /[e (T (p) + e_ (T (p')- 
£m(p+p')]- The coupled integral equations for the renor- 
malized energies must be solved selfconsistently, since 
£cr(p) and ejw(p) appear on both sides of the equations. 

The open-channel collision integral I a takes the usual 
form for a quantum degenerate Fermi gas [12). but 
with a renormalized momentum-dependent interaction 
strength 0|. The collision integrals X a and Tm arise 
from the resonant state coupling and describe the forma- 
tion and dissociation of molecules. They appear in both 
kinetic equations as reciprocal source/sink terms 



4tt 2 k 2 



h* 



dp2 / dp 3 d(p + P2 - P3) 



X6(e<r(p) + £-<t(P2) - ej\/(p3)) 
{[l-/ CT (p)][l-/- ff (p 2 )]/M(p 3 ) 
/<r(p)/-«r(Pa)[l + /af(pa)]}. 



(4) 



The expression for the collision integral Im in JSJ can 
be obtained as Tm — —T-a and by making the following 
substitutions on the r.h.s. of (@J): {p — > Pi,P3 — > p}. 
From this relation it can be shown that the total atom 
population N to t = J dp J drf tot (p,r,t)/h 3 is conserved, 
where / to t(p,r,<) = / T (p, r, t) + f L (p, r, t) + 2/ M (p, r, t). 

Equations through along with the supplemen- 
tary definitions in the text, compose the nonequilibrium 
kinetic theory for a degenerate Fermi gas with resonant 
interactions. In deriving these equations we have made 



two important approximations about the molecules |14| : 
(i) our theory describes bare molecules since we neglect 
the molecular density renormalization and (ii) the finite 
molecular lifetime (given by the imaginary part of the 
molecule self-energy) is neglected in defining the quasi- 
particle energies e CT and 6m- The physics of the lifetime 
of the molecules (molecular dissociation) is, however, de- 
scribed by the collision integrals. The impact of approxi- 
mation (ii) is that it ultimately leads to the appearance of 
the energy conserving delta functions in the collision in- 
tegrals. We note that near the resonance, the molecules 
are actually dressed by the couplingto the atoms and 
arc known to have a size of a(B)/2 OH El, 

which is 

much larger than the size of a bare molecule. We expect 
our treatment to provide a good description when the 
molecules do not overlap spatially: n[a(B)] 3 <C 1. When 
n[a(B)] 3 > 1, the separation of timescales r sc <C t co \ <C t 
assumption at the heart of our theory breaks down. 

It is straightforward to show that the stationary so- 
lutions of (J2J and © are Fermi and Bose distribu- 
tions: /°(p,r) = {e^Mp.r)-^] + an d /° f (p,r) = 
{efSNtP^-H-i}-!, w hcre (3 = l/k B T. In order for T CT 
and Tm to vanish, the chemical potential of the molecules 
must satisfy \x m = Mt + Mi- The collision integral Im de- 
scribes the process by which the molecular component 
comes into chemical and thermal equilibrium with the 
atoms, resulting in a Bose-Einstein distribution for the 
molecules. The molecular population depends strongly 
on the value and time rate of change of the resonant 
state energy e ros [£?(£)]. In our theory, molecules can only 
be produced when the renormalized resonance energy is 
tuned positive; on the negative detuning side of the res- 
onance (e rcs < 0), energy and momentum conservation 
cannot be satisfied, leading to X a = Im = 0. 

Before giving an explicit application of the theory, it 
is important to comment on other relaxation processes 
that are not treated in our present theory. In Fig. 1 we 
give a summary of elastic and inelastic collisions that will 
determine how the atom- molecule mixture relaxes |10( . 
The check marks indicate the processes described by our 
theory and the symbol 7 represents the relaxation rate for 
a given process. There are two other processes that play 
a role in forming molecules: two atoms form a molecule 
under the influence of a third atom (A + A + A «-> M + A) 
or molecule (A + A + M ^ M+M). As discussed above, 
on the negative side of the resonance jaa — 0; however, 
in general ^aaa and jaam do not vanish. Our kinetic 
equations are therefore a good description when ^aaa 
and Jaam are much smaller than the other relaxation 
rates. We finally note that inelastic collisions, such as 
those listed in Fig. 1, are assumed to be negligible. 

We now derive a simplified model to describe the atom- 
molecule population dynamics. We assume an equal mix- 
ture of spin states /f = /j_ = Ja and take Um(?) = 
2Ua{y)- Our goal is to obtain a coupled set of dy- 
namical equations describing the populations Ni(t) and 
energies Ei(t) of the atom and molecule distributions, 
which are defined by iVj(t) — J dr J dp/j(p, r, t)/h 3 and 
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Elastic collisions 


thermal relaxation 


chemical and thermal relaxation 


>/y A ; A + A <-» A+A 


■J"Iaa '■ A + A^-M 


y M : M+M^M+M 


■{ ■„■■ \ ■ -A A \ ■ ■ M ■ \ 


Yam • A + M ^ A + M 


Yaam ■ A A M ■ M .'./ 


Inelastic collisions 


Yam : M- 


-A^M*+A + 


j'aaa ■ A- 


-A + A^M"+A + 



FIG. 1: Summary of different relaxation processes. The letter 
A signifies an atom in either spin component, M refers to 
the resonant molecule, M* describes a more deeply bound 
molecule, and "+" denotes that the particle carries away an 
energy far in excess of the thermal energy IcbT. The check 
marks designate the processes that are treated in this letter. 



Ei(t) — J dr J dpei(p,r,t)fi(p,r,t)/h 3 . The subscript i 
denotes atoms or molecules i G {A, M}. 

Equations of motion for Ni(t) and Ei(t) are obtained 
by differentiating with respect to time and substituting 
the kinetic equations from @ and J3J . We make the fol- 
lowing three simplifications in order to obtain a model 
description: (a) We neglect the renormalized mean field 
contribution to the energies, and everywhere make the 
substitution ej(p, r, t) — > q ' (p, r). (b) For temperatures 
not too deep into the quantum degenerate regime (e.g. 
T/Tp > 0.5 0]), a classical gas provides a reasonably 
good description. We therefore take (1 ± /j) — > 1 in 
the collision integrals, (c) We finally assume Maxwell- 
Boltzmann forms for the distributions, but allow the 
temperatures and chemical potentials to vary in time: 
f A (p,r,t) = exp{-[e^(p,r) - /j, A (t)]/k B T A (t)} and 

/ M (p,r,i) = exp{-[e^ } (p,r) - ^M(t)}/k B T M (t)}. 

It is now straightforward to obtain a closed set of cou- 
pled differential equations for the populations Ni(t) and 
energies Ei(t). The atom equations are 



dN A 

dt 
dE A 

dt 



Treves) 

7«(erc S ) 



N 



M 



N A e 



-(t Ias -fl,A)/k-BTA 



(5) 



E 



M 



(E A + e ICS N A )e-^-^/ k * T * 



(6) 



The variables (ii(t) and Ti(t) are related to the system 
variables {N A (t),N M (t),E A (t),E M (t}} using the defi- 
nitions of Ni(t) and Ei(t) given above. Due to num- 
ber conservation, the molecule population is given by 
-/Vm(£) = Nt t/2 — N A (t). The equation for the molecule 
energy is dEM(t)/dt = e Tes NM(t) — 2dE A (t)/dt, where 
e rcs is the time rate of change of the resonant state en- 
ergy. The total energy E tot (t) = 2E A (t) + E M {t) varies 
as dEtot/dt = e res NM(t). It can be shown that the r.h.s. 
of JSJ) and © vanish when (j,m = 2\i A and Tm = T A . 
The relaxation rate 7 K (e) is given by 



7«W 



27r?i 4 



(7) 



which has the form of a Fermi's golden rule |?J. The 
Heaviside step function 9(e) signifies that for negative 
detunings, the exchange of atoms and molecules ceases. 
This rate determines whether a sweep of e ros (t) is adia- 
batic in the sense that the entropy is conserved. Using di- 
mensional analysis, we find a good measure of adiabatic- 
ity in our calculations is a/e rcs , where a = k B T ~f K (k B T). 
The process is adiabatic if a/i rcs 3> 1. 

To express this inequality in terms of physical param- 
eters, we approximate the coupling strength by k 2 w 
27r?i 2 abgA/i AB/m 16.], where a^ g is the background 
scattering length, Afj, is the difference in magnetic mo- 
ments of the open and closed channels, and AB is the 
resonance width. Our theory does not provide the mag- 
netic field dependence of e Ies (B), so we make the sim- 
ple approximation e res ~ A/iB. With these estimates, 
we obtain a/e res w (2Tr) 3 ^ 2 (hab g /m\^ h )AB/B, where 
Ath = h/ \J2 r Kmk B T is the thermal deBroglie wavelength. 
This expression resembles the Landau-Zener adiabaticity 
requirement jfj], but with the volume factor A^ h appear- 
ing for a thcrmalized gas at finite temperatures. 

We solve the coupled population and energy equations 
numerically using physical parameters for the recent ex- 
periments by Regal et al. 0, Q with 40 K at the reso- 
nance B = 20.21 mT, for which AB = 0.78 mT and 
a bg = 174ao, where a is the Bohr radius. We take 
Ntot = 10 6 and trap frequencies v x = v y = 290 Hz and 
v z = v x /<i0- The Fermi temperature is Tp = 0.47//K. 
We present results for T/Tp = 0.5, which is in the range 
0.1 < T/Tp < 2.0 that applies to recent experiments on 
molecule formation Taking A^i w 2^b, the rate con- 
stant is ct/fiB — 1-9 T/s, where /ib is the Bohr magneton. 

The resonance energy is ramped down linearly as 
e res (i) = e r cs(£ = 0) — e rcs t and stopped on the negative 
detuning side. In Fig. 2, we plot the molecule fraction 
t]m = 2NM/N to t versus e xes (t) for two different sweep 
rates: a/e ros = 5 (red dashed line) and a/e rcs = 0.2 
(blue dot-dashed line). For comparison, we plot the equi- 
librium solution for a classical ideal gas mixture of atoms 
and molecules along a line of constant entropy, given by 
the solid black line. When the adiabaticity criterion is 
satisfied a/e res 3> 1, the sweep achieves the maximum 
transfer allowed by following along a line of constant en- 
tropy, before entering the negative side of the resonance 
where the atom-molecule transfer stops. Recall our as- 
sumption that the rates "f AAA and 7aam appearing in 
Fig. 1 are much smaller than j AA ~ j K . If the ramp is 
carried out slowly compared to these 3-body relaxation 
rates, the system should follow the solid black curve on 
the negative side, which approaches unity. 

In Fig. 3 we plot the final atom fraction r] A = \ — t\m 
versus the inverse sweep rate. An important result is that 
it is impossible to achieve a perfect transfer of all the 
atoms into molecules at finite temperatures, even when 
the sweep is adiabatic; one can only achieve the limiting 
value minl?^}. This saturation of the transfer efficiency 
is observed in experiments 0,0- The value of minl?^} is 
strongly dependent on temperature, as seen in the inset of 
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FIG. 2: Molecule fraction versus the resonance energy for two 
different sweep rates. The red dashed and blue dot-dashed 
curves are obtained by integrating the equations of motion for 
{NA(t),N M (t),EA{T),E M (t)} with a linear ramp of e res (t), 
starting with e las (t — 0)//cb = 2 /iK. The black solid curve 
is obtained by calculating the equilibrium solution for a clas- 
sical ideal gas mixture of atoms and molecules. The initial 
temperature is T(t = 0)/T F = 0.5. 



Fig. 3 where we plot minf?^} versus T. Our model pre- 
dicts that the transfer efficiency (1 — minl?^}) increases 
as T decreases. We note that for temperatures below 
T/T-p ps 0.5, the effects of quantum statistics, which are 
neglected in our simple model, will become important. 

We have presented kinetic equations and de- 
scribing a resonantly coupled atom-molecule mixture in 
the pseudogap regime above the superfluid transition. As 
a first application, we derived a non-interacting model for 
the population dynamics, similar to the toy model pre- 
sented in Ref. |l9j for equilibrium properties. The satu- 
ration of molecule formation predicted by our model was 
observed in recent experiments above yj] and below 
the superfluid transition temperature. We also predict 
that the transfer efficiency increases as the temperature 
is lowered. The kinetic equations provide a foundation 
for many future applications of dynamical effects, such 
as exploring the collisionlcss to hydrodynamic crossover 
of the damped collective excitations. 

We acknowledge J. Cooper, A. Griffin, M. Holland, P. 
Julienne, and E. Tiesinga for helpful discussions. 
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FIG. 3: Atom fraction versus the inverse sweep rate. The 
blue circle and red cross correspond to the blue dot-dashed 
and red dashed curves shown in Fig. 2, respectively. The inset 
shows how min{77A} decreases with decreasing temperature. 
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